IM>-f(lS7  842 
UNCLASSIFIED 


STUDV  ON  RLGORITHHS  FOR  PREDICTION  OF  SOLID  PROPELLANT 
ROCKET  NOTOR  PERFORHRNCE(U>  FOREIGN  TECHNOLOGY  DIV 
HRIGHT-PATTERSON  AF8  OH  X  HANG  ET  AL.  22  APR  8G 
FTD-I0(RS>T-0192-8£  F/0  21/8.  2 


1/1 


AD"A167  842 


I  “  .  -  I  -  —  — 

FTD-ID(RS)T-0152-86 


FOREIGN  TECHNOLOGY  DIVISION 


STUDY  ON  ALGORITHMS  FOR  PREDICTION  OF  SOLID  PROPELLANT  ROCKET  MOTOR 

PERFORMANCE 

by 

Wang  Xln  and  Zhang  Zhongqln 


Approved  for  public  release; 
Distribution  unlimited. 


PTO-  lD(RS)T-0152-86 


HUMAN  TRANSLATION 

FTD-ID(RS)T-0152-86  22  April  1986 

MICROFICHE  NR:  PTIV86-C-001752 

STUDY  ON  ALGORITHMS  FOR  PREDICTION  OF  SOLID  PROPELLANT 
ROCKET  MOTOR  PERFORMANCE 


By:  Hang  Xin  and  Zhang  Zhongqin 
English  pages:  12 

Source:  Hang  Kong  Xuebao«  Vol  5,  Nr.  3,  September  1984, 
pp.  288-295 

Country  of  origin:  China 
Translated  by:  FLS,  INC. 

F33657-85-D-2079 
Requester:  FTD/SDBP 

Approved  for  public  release;  Distribution  unlimited. 


THIS  TRANSLATION  IS  A  RENDITION  OF  THE  ORIGI- 

NAL  FOREIGN  TEXT  WITHOUT  ANY  ANALYTICAL  OR 

PREPARED  BY; 

EDITORIAL  COMMENT.  STATEMENTS  OR  THEORIES 

ADVOCATED  OR  IMPLIED  ARE  THOSE  OF  THE  SOURCE 

TRANSLATION  DIVISION 

AND  DO  NOT  NECESSARILY  REFLECT  THE  POSITION 

FOREIGN  TECHNOLOGY  DIVISION 

OR  OPINION  OF  THE  FOREIGN  TECHNOLOGY  DIVISION 

•  WPAFB,  OHIO. 

FTD-  ID(RS)T-0152-86 


Pate  22  April _ 19  86 


GRAPHICS  DISCLAIMER 

All  figures,  graphics,  tables,  equations,  etc.  merged  into  this 
translation  were  extracted  from  the  best  quality  copy  available. 


I  Accession  Foi^ 

'l.'TIS 

nr  1C  T/.H 

U;'-'iunoi!!iced 

Ja::t  if  icntioiu 


□ 


Fy- 


'  /vni  1 1 Codes 
1  ij'vr  U.  and/or 


^ ,c lal 


STUDY  ON  ALGORITHMS  FOR  PREDICTION  OF  SOLID  PROPELLANT  ROCKET  MOTOR 
PERFORMANCE 


Wang  Xin  and  Zhang  Zhongqin 


Bei jing 
Recei ved 


Institute  o-f  Aeronautics  and  Astronautics 
on  ^eot  ember  17,  1986 


This  article  introduces  two  predictive 
algori th%^-for  the  per-formance  o-f  solid  propellant 
rocket  motor.  The  emphasis  is  on  the  introduction 
of  Time-Space  Algorithm.  The  authors  proposed  a 
general  two-dimensional  grain  calculation 
procedure  in  order  to  conduct  the  grain 
calculations.  Therefore,  the  predictive  algorithms 
for  performance  introduced  by  this  article  show 
generality.  A  comprehensive  computer  program  for 
the  aforementioned  method  has  been  written  and 
applied  to  calculating  the  performance  of  three 
different  solid  rocket  motors.  The  calculated 
results  are  consistent  with  those  derived  from 
experimental  data.^ 


Nomencl ature 


a  speed  of  sound 

combust!  on  speed  coefficient 
Aacritical  speed  of  sound 
combustion  area 
^combustion  area  at  the  head 
passage  area 
A,  nozzle  throat  area 
C» exhaust  gas  coefficient 
h  space  increment 
H enthalpy 

^^enthalpy  of  propellant 
m  mass  flow  rate 
•  pressure  index 
P  combustion  gas  pressure 


f combustion  speed  of  propellant 
K gas  constant 

7*  combustion  gas  temperature 
9  circumference  of  combustion 
T*adiabatic  combustion  temperature 
I  time 

*  combustion  gas  speed 
X  axial  coordinate 
y  specific  heat  ratio 
A  speed  coefficient 
P  combustion  gas  density 
Prpropellant  density 
«  ratio  of  corrosion 
T  time  increment 


Subscript 


0  parameters  at  0  cross-section 
of  grain 

/  parameters  at  passage  exit 
cross-sect i on 


I  parameters  at  nozzle  throat 
$  stagnation  parameter 
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I.  Preface 

The  purpose  for  predicting  the  performance  of  a  solid  propellant 
rocket  motor  is:  to  base  on  raw  data  such  as  structure  of  the  motor, 
geometrical  shape  and  dimensions  of  the  grain,  characteristics  of 
propel 1  ant , etc.  and  calculate  the  changes  with  respect  to  time  and 
their  distribution  along  the  length  of  the  combustion  chamber  for 
parameters  such  as  combustion  gas  pressure  flow  speed  9,  density P 

and  temperature  T  during  motor  operating  period. 

Great  advances  have  been  achieved  in  the  development  of  performance 

algorithms  throughout  the  years.  In  earlier  days  the  performance 

algorithms  were  built  on  the  basis  of  quasi-steady  flow  and 

zero-dimension  simplification  Cil.  Then,  a  one-dimensional  model  was 

adopted  for  conditions  with  larger  volume  filling  coefficient  and 

rather  distinctive  corrosion  combustion  phenomenon.  A  group  of 

ordinary  differential  equations  C2,51  could  be  obtained  when 

conducting  quasi-steady  simplification  on  the  basis  of  a  one-dimensional 

model,  and  generally  they  could  be  solved  using  the  Runge-Kutta  or  the 

Merson  method  C73.  If  the  unsteady  effects  are  considered  on  the  basis 
the 

of Aone-dimensional  model ,  then  the  group  of  performance  control 
equations  are  of  first-order,  linear  hyperbolic  partial  differential 
equations,  and  the  finite-difference  and  the  characteristic-line 
method  could  be  applied  to  solving  them  [83. 

This  article  concentrates  on  the  introduction  of  general  algorithms 
for  quasi-steady  and  unsteady  equations,  which  is  the  Time-Space 
Alternate  Algorithm.  It  was  developed  on  the  basis  of  the  Triple-Loop 
Iteration  Algorithm.  Since  the  grain  and  performance  calculations  are 
closely  related  and  indivisible,  a  general  two-dimensional  grain 
calculation  procedure  was  employed  in  these  calculations. 

II.  Triple-Loop  Iteration  and  Time-Space  Alternate  Algorithms 

1.  Motor  Model  and  Its  One-Dimensional  Performance  Control  Equations 


Select  motor  model  as  depicted  in  Figure  1 


The  -following  assumptions  are  made  -for  the  -flow  field  inside  the 
combustion  chamberi  combustion  gas  is  ideal  gas;  flow  is 
one-dimensional  in  the  grain  passage; flow  process  is  frictionless, 
adiabatic,  chemically  nonreactive  and  homogeneous  gaseous  flow;  the 
x-axis  component  of  momentum  of  combustion  gases  generated  on  the 
combustion  surface  is  zero.  The  performance  control  equations  can  be 
derived  as: 
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If  :•«<»,  P<CP»,  (within  combustion  gases  resident  time),  the 
operation  process  of  motor  can  be  considered  quasi -steady  and  the 
above  equations  can  be  simplified  as: 
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Figure  1.  Schematic  of  solid 
propellant  rocket  motor 


Iteration  Algorithm 

flow,  the  energy  equation  (7)  in  the  one-dimensional 

equations  (5) -(6)  can  be  simplified 

(») 


Then,  only  the  continuity  and  momentum  equations  in  the  group  are  in 
differential  form.  Now  let's  select  the  mass  flow  rate  m  and  P pressure 
as  var i ables «and  change  the  continuity  and  momentum  equations  into  the 


-following  calculation  formal 
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Every  unknown  parameter  can  be  expressed  as  a  function  of  m  AntiP 
using  the  energy  and  state  equations  and  other  well  known  formulae. 


The  solving  of  performance  equations  is  a  boundary  value  problem 
rather  than  an  initial  value  one;  therefore  the  test-firing  method  is 
adopted.  First,  assume  that  the  test  value  of  head  presure  Pt  can  be 
estimated  according  to  the  following  formula: 


where 
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and  Xi  is  determined  by  the  following  formula* 

The  test  value  of  P»  at  time  #,  can  be  obtained  by  selecting  the  A 
value  that  satisfies  the  flow  conditions  at  time  t,  .  As  to  the  Pt  test 
values  at  time  f»  and  beyond, they  can  be  obtained  by  using  the  P»  values 
at  two  previous  time  instants  and  applying  the  proportionality  method, 
which  is: 

P.(/.)-P.«-..)-^{^;5-.  (»»-2.  3.  •)  (15) 


At  any  time  instant,  the  mass  flow  rate  across  the  0 
should  be  equal  to  the  combustion  gas  generation  rate 
head,  then 


cross-section 
at  the  motor 

(16) 


Once  the  grain  head  boundary  conditions  Pt  and  are  determined, 
equations  (10)  and  (11)  can  be  solved  as  initial  value  problems  using 
the  Merson  method.  Then  adjust  Pt  value  according  to  the  boundary 
conditions  at  the  grain  tail.  In  order  to  do  this,  the  combustion  gas 


generation  rate  <«»  and  exhaust  rate  must  be  obtained. 


1X4  • 

m, 


Mhere  is  the  combustion  gas  generation  rate  at  the  cross-section  o-f 
grain  tail.  Since  the  -flow  -from  the  cross-section  o-f  grain  exit  to 
nozzle  throat  is  also  considered  quasi -steady  •«,  should  be  equal  to  . 
If  the  calculated  values  for  and  are  not  equal,  then  this  shows 
that  the  assumed  value  for  P,  is  not  reasonable  and  should  be  adjusted. 
The  proportionality  method  can  be  used  for  the  first  adjustment  of  the 
P»  test  value,  that  is: 


(17) 


After  the  number  of  iterations  exceeds  two,  interpolation  formula  can 
be  used  to  adjust  Po  value.  At  this  time,  let  A»»« then 
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The  iteration  continues  until  A'n  is  less  than  a  designated  value. 

When  the  solution  at  time  has  been  obtained  and  the  solution  at 
time  tm  is  continued,  the  geometrical  parameters  at  time  tm  should  be 
obtained  by  grain  calculations  first.  The  grain  thickness  burned  at 
each  node  during  time  At  can  be  obtained  according  to  the  following 
calculation: 

where  n'is  the  average  speed  of  combustion  during  time  At  and  it  is 
dependent  upon  the  speed  of  combustion  at  time  ^  .  However,  the  speed 
of  combustion  at  time  I.  is  affected  by  combustion  gas  parameters  and 
thus  by  Ayi ,  it  must  be  calculated  byjf*^teration  method.  We  employed  the 
scheme  of  simple  prediction  and  multiple  adjustments.  The  specific 
formulae  are: 

rr(W-h«-.i)»  ("-1.  t*  •••) 
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In  applying  the  Merson  method  to  conduct  numerical  integration,  the 
grain  geometrical  parameters  and  their  derivatives  with  respect  to  x 
at  certain  points  between  two  nodes  must  be  used.  These  data  are 
obtained  from  the  LaGrange  quadratic  polynomial  interpolation  and 


their  corresponding  numerical  dif -ferential  equations. 


3.  Time-Space  Alternate  Algorithm 

The  one-dimensional,  unsteady  per-formance  control  equations  (l)-<4) 
are  o-f  first-degree,  near-linear,  hyperbolic  partial  differential 
equations.  If  the  characteristic-line  method  is  used  for  solution,  the 
time  increment  must  be  set  very  small  due  to  the  limitation  of 
stability  conditions  of  the  difference  scheme.  Ten  thousand  or  more 
calculation  steps  are  required  for  a  typical  motor.  This  is  very 
inconvenient  and  sometimes  even  unacceptable. 


The  essence  of  the  Time-Space  Alternate  Algorithm  is  to  reduce 
partial  differential  equation  to  ordinary  differential  equation.  High 
accuracy  for  time  and  space  increments  can  be  obtained.  In  reference 
C61,Burstein  used  a  similar  method  to  solve  hyperbolic  equations  and 
obtained  third-degree  accuracy  for  both  time  and  space  increments. 
This  article  applied  the  Time-Space  Alternate  Algorithm  and  obtained 
not  only  fourth-degree  accuracy  for  time  and  space  increments,  but 
also  nicely  combined  the  solving  methods  of  unsteady  and  quasi -steady 
equations. 


In  the  Time-Space  Alternate  Algorithm,  p,  » and  P  were  selected  as 
caculation  variables.  First  of  al 1 , equations  (l)-(3)  were  changed  to; 


where 
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Then  equations  (19) -(21)  were  written 
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calculation  forms, 


(22) 


thus 


(24) 


at  any  prescribed  time  instant,  i-f  the  values  of  4P/M,  and  *p/M  in 
equations  (22) -(24)  have  already  been  obtained.  Thus,  the  equations 
can  be  considered  ordinary  differential  equations  with  respect  to  x. 

With  this  in  mind,  the  solutions  can  be  obtained  by  the  aforementioned 
Triple-Loop  Iteration  Algorithm.  After  obtaining  the  values  of  p,  o, 
and  #  in  the  axial  direction  at  each  node  and,  respectively, 
substracting  from  them  the  values  of  P  ,  *  and  P  at  the  same  node  but 
at  a  previous  time  instant,  and  then  divided  by  the  time  increment 
between  the  two  time  instants,  the  average  rate  of  change  of  p ,  e  and P 
at  any  node  between  two  time  instants  can  be  obtained.  Using  these 
average  rates  of  change  as  the  approximate  values  for  fP/f#,  andfp/^ 

and  substituting  them  into  equations  (22)-(24)  for  solution  based  on 
the  above  algorithm.  Thus,  repetitive  and  alternating  calculations  can 
be  conducted  according  to  time  and  space  increments  so  as  to  adjust 
the  values  of  fP/ff,  fv/M ,  and  fp/^  i  terati  vel  y .  Once  a  certain  accuracy  is 
satisfied,  the  calculations  for  the  next  time  instant  begin  and 
proceed  till  the  motor  stops. 


If  equations  (22)- (24)  are  expressed  in  vector  forms,  then 


A 


C25) 


When  discussing  the  stability  of  the  difference  scheme,  the  above 
equation  is  usually  simplified  as: 


(26) 


where  is  a  constant  matrix. 

When  using  the  Time-Space  Alternate  Algorithm  to  solve  model 
equation  (26) ,  it  is  equivalent  to  applying  the  difference  scheme: 


Here  the  di-f-ference  scheme  is  unconditionally  stable  with  respect  to 
the  model  equation.  Through  actual  calculations,  it  is  proven  that  -for 
actual  equations,  the  above  discussions  also  apply. 

The  Triple-Loop  Iteration  Algorithm  was  adopted  in  consideration  o-f 
the  solving  o-f  equations  with  respect  to  x  in  the  Time-Space  Alternate 
Algorithm.  There-fore,  only  the  portion  di-f-ferent  -from  the  algorithm 
-for  quasi-steady  equations  is  introduced  here. 

The  necessary  initial  conditions  -for  the  Time-Space  Alternate 
Algorithm  include  the  values  o-f  P,  »  and  P  at  each  node  at  the  initial 
time  instant  and  the  values  o-f  4P/m,  »pfM  and  *olM.  When  conducting 
calculations  immediately  -following  the  startup  o-f  the  motor,  the  above 
initial  conditions  could  be  determined  by  the  calculations  -for  the 
startup  process.  However,  when  the  calculations  were  started  directly 

at  the  steady  state  o-f  the  motor,  we  assumed  that  the  -flow  -field 

a 

inside  the  combustion  chamber  reached  steady  state  st  the  initial 
instant;  thus  the  values  o-f  *P/^,  and  fpM  at  each  node  can  be 

assigned  zero  and  the  values  for  P,  •  and  p  at  each  node  obtained  by 
the  Triple-Loop  Iteration  Algorithm  were  used  as  the  initial  values. 

The  determination  of  the  value  P%  in  the  grain  head  boundary 
conditions  is  the  same  as  that  of  the  Triple-Loop  Iteration  Algorithm. 
But  the  value  of  P#  and  is  determined  by  using  the  following  formula: 

P,«p,//?r » 

At  time  ,  the  test  values  for  *i>/M  And  *p/*t  values  at  each 

node  are  substituted  by  the  average  values  AP/Af,  Av/A(  and  Ap/A-’  between 
time  «  and  f—i  .  After  the  values  of  P  ,  v  and  P  at  each  node  at  timef 
were  obtained  they  are  then  adjusted.  The  adjustment  formulae  are: 

\  ^  - 

(m-J,  2,  — »  B-J,  2,  — ) 

Adjustment  can  be  conducted  repetitively  until  the  difference  of  P ,p 
and  P  values  between  two  consecutive  iterations  is  smaller  than  a 
prescribed  value. 


III.  Grain  Calculation 


The  purpose  -for  grain  calculation  is  to  obtain  the  grain 
geometrical  parameters  required  -for  the  motor  performance  calculation. 
For  different  types  of  motor  performance  calculation,  not  only  the 
required  grain  geometrical  parameters  are  different,  the  type  of 
coordination  between  the  grain  and  motor  performance  calculation  is 
also  different.  For  zero-dimensional  performance  algorithm,  the 
required  grain  geometrical  parameters  generally  correspond  to  the 
value  of  grain  combustion  surface  of  a  given  burned  grain  thickness. 
Then,  the  grain  and  performance  calculations  can  be  conducted 
separately.  However,  for  the  performance  algorithm  of  one-dimensional 
and  variable  cross-sectional  passage  the  required  grain  geometrical 
parameters  generally  correspond  to  the  values  of  the  combustion 
circumference  of  the  local  burned  grain  thickness  and  passage  surface 
area.  With  corrosion  combustion  present,  the  grain  and  performance 
calculations  are  already  closely  related  and  indivisible,  and  they 
must  be  conducted  alternately  through  iteration. 

The  performance  algorithm  introduced  in  this  article  incorporates 
the  grain  calculations.  For  this  reason  the  authors  proposed  a 
two-dimensional  general  algorithm  for  grain  calculation  C93.  To  expand 
from  the  basis  of  this  algorithm,  three-dimensional  calculation 
problems  such  as  the  double-arc  transient  head  seal  and  the  tail-end 
surface,  which  is  located  at  the  ellipsoid  behind  the  motor  and  does 
not  restrict  combustion. 

IV.  Examples  and  Discussions 

Test  runs  of  this  computer  program  on  three  different  real  motors 
were  conducted.  Both  Triple-Loop  Iteration  and  Time-Space  Alternate 
Algorithms  were  applied  on  each  motor.  For  example  1  and  3,  since 
their  performance  parameters  during  the  operation  period  changed 
little  with  respect  to  time  and  the  results  obtained  from  both 
algorithms  differed  little,  only  the  results  obtained  from  the 
Triple-Loop  Iteration  Algorithm  are  presented.  As  for  Example  2,  the 
results  obtained  from  both  algorithms  are  presented. 


Example  1.  The  grain  is  star — shaped,  two-dimensional  and  internal 
cavi ty— f i 1 led.  The  head  portion  is  double-arc  transient  head  seal.  The 
tail-end  surface  is  located  at  the  ellipsoid  behind  the  motor  and 
allows  unrestricted  combustion.  See  Figure  2  for  the  results  of 
calculation. 

Example  2.  The  grain  is  wheel-shaped  and  internal  cavity-filled. 
Both  the  head  and  tail-end  surfaces  are  located  at  the  cylindrical 
section  of  the  motor  and  both  allow  unrestricted  combustion.  See 
Figure  3  for  the  results  of  calculation. 

Example  3.  The  grain  is  of  inner/outer  surfaces  simultaneous 
combustion  and  single  tubular  fill.  Both  the  head  and  tail-end 
surfaces  restrict  combustion.  See  Figure  4  for  results  of  calculation. 


calculated  curve  (Triple-Loop 


Figure  2.  Pressure  vs  time  diagram  of  solid 
rocket  motor  for  Example  1. 

It  can  be  observed  from  the  calculated  results  that  when  changes 
with  respect  to  time  in  the  performance  parameters  are  more 
significant  (Example  2) ,  there  are  marked  differences  between  results 
obtained  from  the  quasi-steady  ordinary  differential  equations 
(Triple-Loop  Iteration  Algorithm)  and  those  from  the  unsteady  ordinary 
differential  equations  (Time-Space  Alternate  Algorithm) ,  and  the 
latter  had  a  calculation  accuracy  about  1.67.  higher  than  that  of  the 
former . 
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Figure  3.  Pressure  vs  time  diagram  of  solid 
rocket  motor  for  Example  2. 
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Figure  4.  Pressure  vs  time  diagram  of  solid 
rocket  motor  for  Example  3. 
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STUDY  ON  ALGORITHMS  FOR  PREDICTION  OF 
SOLID  PROPELLANT  ROCKET  MOTOR  PERFORMANCE 

Wang  Xbt  and  Zhang  Zhongqin 
iBeijing  Insiitidt  of  AeronaUics  and  Astronautics) 

Abstract 

Two  algorithms  for  prediction  of  solid  propellant  rocket  motor  perfor¬ 
mance  are  introduced.  Triple  Loop  Iteration  Algorithm  it  applied  to  solu¬ 
tion  of  quasi-steady  equationsC  5  )‘^(  8  ).  Its  computer  program  comprises 
a  procedure  with  three  iterative  loops.  Time-Space  Alternate  Algorithm  can 
be  applied  to  soluting  both  quasi-steady  equations  and  unsteady  equations 
(22)'^(24).  If  »v/at,  ap/it  in  equations  (22)~(24)  are  gixen,  equa¬ 

tions  (22)^(2t)  will  be  reduced  to  ordinary  differential  equations.  They 
can  be  solved  by  Merson's  method.  When  AP/At,  Av/At,  Af>/At  are  appro¬ 
ximated  to  iP/at,  a-jjat,  apfat,  this  algorithm  for  model  equation 

eg  . 

as 

corresponds  to  the  finite  difference  scheme  v 


A  f  aJ"  ij*  »r*  'i 


It  is  an  unconditionally  stable  scheme. 

A  general  two-dimensional  grain  calculation  method  is  introduced  for 
calculation  of  grain.  A  comprehensive  computer  program  for  above  men¬ 
tioned  methods  has  been  written  and  applied  to  calculating  the  perfor¬ 
mance  of  three  different  motors.  The  results  of  the  numerical  examples 
a.e  consistent  with  the  experimental  data. 


